Electron-phonon coupling in quasi free-standing graphene 
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Quasi free-standing monolayer graphene can be produced by intercalating species like oxygen or 
hydrogen between epitaxial graphene and the substrate crystal. If the graphene is indeed decoupled 
from the substrate, one would expect the observation of a similar electronic dispersion and many- 
body effects, irrespective of the substrate and the material used to achieve the decoupling. Here 
we investigate the electron-phonon coupling in two different types of quasi free-standing monolayer 
graphene: decoupled from SiC via hydrogen intercalation and decoupled from Ir via oxygen interca- 
lation. Both systems show a similar overall behaviour of the self-energy and a weak renormalization 
of the bands near the Fermi energy. The electron-phonon coupling is found to be sufficiently weak 
to make the precise determination of the coupling constant A through renormalization difficult. The 
estimated value of A is 0.05(3) for both systems. 



The unusual electronic structure of graphene (TJ [2] has 
significant implications on the transport properties as 
well as on the many-body effects in this material. An 
interaction of particular importance for applications is 
the electron-phonon coupling, as this mechanism can be 
dominant for the scattering of carriers near the Fermi 
surface. In theoretical work, it was recognised early that 
the Dirac-like dispersion and its corresponding restriction 
of scattering phase space renders the electron-phonon 
coupling in graphene quite different from that in ordi- 
nary metals O |4]. In particular, it leads to a doping- 
dependent electron-phonon coupling strength A and a 
much more complex shape of the self-energy than in or- 
dinary metals, where the density of states near the Fermi 
level, and its contribution to the self-energy, is approxi- 
mately constant on a phonon energy scale. 

Angle-resolved photoemission (ARPES) is a particu- 
larly valuable technique to address many-body interac- 
tions experimentally since it provides detailed informa- 
tion on the energy and /c-dependent lifetime of the carri- 
ers. Several ARPES studies of the electron-phonon cou- 
pling in graphene have been published [5HT2] , mostly for 
strongly n-doped graphene (n ~ 10''^^ cm~^), where the 
electron-phonon scattering was found to be of intermedi- 
ate strength with A in the order of 0.2 — 0.3. 

For weakly doped graphene, the reduction of scatter- 
ing phase space implies a weaker and eventually vanishing 
electron-phonon coupling strength. Experimentally, it is 
difficult to verify this with high precision because of the 
often employed methodology to determine A. This relies 
on the observed energy dependence of the electronic self- 
energy near the Fermi energy Ep^ where electron-phonon 
coupling gives rise to kinks in the dispersion and a de- 
crease of the linewidth. For very weak coupling, these 
spectral features are equally weak and hard to distin- 
guish from the noise in the spectra. It should also be 
mentioned that the determination of A in this way only 



works for sample temperatures that are low compared 
to the typical excitation temperature for phonons. This 
is seemingly easy to assure for graphene with its high 
Debye temperature, but it has been questioned if the 
Debye temperature or the doping-dependent and lower 
Bloch-Griineisen temperature sets the relevant scale for 
electron-phonon coupling in graphene [13 . Neverthe- 
less, there is mounting experimental evidence that the 
electron-phonon coupling in weakly doped graphene is 
indeed quite small [[TTl [Til [T5] . 

An essential question for any meaningful comparison 
between experimentally determined many-body effects 
and calculations is of course how much the properties 
of the epitaxial graphene, typically used in experiments, 
actually resemble those of the pristine graphene assumed 
in calculations. Evidence from transport measurements 
suggest that the contact of graphene to a substrate can 
significantly reduce the carrier mobility, especially at 
higher temperatures [lH [17]. Epitaxial graphene on 
both SiC and transition metal surfaces still exhibits a 
relatively strong interaction with the substrate, at least 
strong enough for the graphene to have a well-defined az- 
imuthal orientation with respect to the substrate. In rare 
cases, such as for Pt (111) [18], this well-defined azimuthal 
orientation is lost, but for the majority of samples even 
stronger interactions are directly observed in the disper- 
sion of the graphene bands, ranging from weak replica 
bands with mini-gaps on Ir(lll) [19^,^ to interactions 
which are so strong that the typical graphene dispersion 
is not observed, even though the structure of the layer 
is graphene-like [27, "2? . Apart from these largely struc- 
tural effects, a relatively strong doping is also frequently 
observed. Note that even in the absence of a graphene- 
substrate interaction immediately evident in the band 
dispersion, the different dielectric properties of the sub- 
strate (wide-gap semiconductor vs. metal) might still 
affect the many-body interactions in graphene. 
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In this paper, we investigate the so-cahed quasi free- 
standing monolayer graphene (QFMLG) as a possible 
experimental solution to the problem of the graphene- 
substrate interaction. QFMLG was first proposed in con- 
nection with hydrogen-intercalation below the graphene- 
like so-called buffer layer forming the (6a/3 x 6a/3)R30° 
structure on SiC(OOOl) 03]. Whereas the buffer layer 
does not show a graphene-like electronic structure [21 , a 
weakly p-doped and narrow Dirac cone emerges upon hy- 
drogen intercalation [24 . Similar results can be achieved 
by oxygen intercalation for graphene on Ir(lll) [25] 
where the graphene is also p-doped but any other in- 
dication of graphene-substrate interaction is lost. Here 
we probe the electron-phonon coupling on these two 
QFMLG systems that have been synthesised via hydro- 
gen intercalation on SiC(OOOl) and oxygen-intercalation 
on Ir(lll). We find that the electron-phonon coupling is 
weak A ^ 0.07 for both systems. In view of the rather 
different substrates and materials used for intercalation, 
this is taken as a suggestion that the QFMLG can indeed 
be taken as a model for pristine graphene. 

High-quality graphene was prepared using standard 
methods on Ir(lll) [26j. The sample quality was checked 
with low-energy electron diffraction, showing an intense 
moire pattern for a clean graphene monolayer. Inter- 
calation of oxygen was achieved by placing the sample 
within a custom made high-pressure O2 doser and main- 
taining a background O2 pressure at 4 x 10~^ mbar for 10 
min while keeping the sample temperature at 523 K. H- 
intercalated QFMLG on SiC(OOOl) was prepared ex-situ 
by the methodology outlined in Ref. In the ARPES 
vacuum chamber, these latter samples were cleaned by 
annealing to 673 K. The ARPES measurements were car- 
ried out at the SGM-3 beamline of the synchrotron radi- 
ation source ASTRID [28 under ultrahigh vacuum con- 
ditions with a base pressure in the 10~^^ mbar range and 
with the sample temperature kept at 70 K. The employed 
photon energies were 47 eV and 32 eV for graphene on Ir 
and SiC, respectively, and the total energy and k resolu- 
tion were better than 18 meV and 0.01 A~^, respectively. 

The measured spectral function for QFMLG on both 
substrates is shown in Figure [l] The data are cuts 
through a three-dimensional data set (taken as a function 
of electron kinetic energy and two emission angles) close 
to the Fermi energy Ep and near the ^-point of the Bril- 
louin zone. Note that here the K-point is defining the ori- 
gin of the /c-coordinate system. The characteristic Dirac 
cone is easily identified. The observed electronic struc- 
ture is drastically different from the situation without 
intercalated oxygen or hydrogen: On Ir(lll) the replica 
bands are not observed any more and on SiC (0001) the 
intercalation is required for any observable Dirac cone. In 
both cases, we find a clear p doping with estimated carrier 
densities of 4 x 10^^ cm~^ (Ir) and 5 x 10^^ cm~^ (SiC). 
On Ir(lll), the p doping is ascribed to a charge transfer 
to the intercalated oxygen. On SiC (0001), it is caused by 
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FIG. 1: Electronic structure of QFMLG on Ir(lll) (a,c) and 
SiC(OOOl) (b,d) measured by ARPES. Shown is the photoe- 
mission intensity as a function of k and binding energy (on 
the top part) and as a function of both components of k at 
the Fermi energy (on the bottom part). Dark corresponds to 
high intensity. The insets show the direction of the shown 
dispersion within the Brillouin zone. 



a spontaneous polarization of the substrate [29^. 

Data as in Figure [l] form the basis for a more detailed 
analysis of the electron-phonon coupling. For QFMLG 
on Ir(lll), this analysis is carried out using the disper- 
sion perpendicular to the T — K direction. For QFMLG 
on SiC, the chosen direction is not a high-symmetry di- 
rection but encloses an angle of 40.4° with the T — K 
direction. Carrying out the same analysis for other di- 
rections in /c-space does not lead to a significantly dif- 
ferent result, consistent with the theoretical expectation 
of a uniform electron-phonon coupling strength over the 
Fermi contour [8 . 

Information about the electron-phonon coupling needs 
to be extracted in two steps. The first is to obtain the 
self-energy from the data and the second to extract infor- 
mation about the many-body effects from the self-energy. 
Assuming no /c-dependence of the self-energy S and no 
strong matrix element effects, ARPES measures the spec- 
tral function 
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where huo is the binding energy, e{k) the bare dispersion 
(without many-body effects) and S the self-energy with 
real part and imaginary part For a cut at constant 
energy, a so-called momentum distribution curve (MDC), 
this has the form of a Lorentzian line with the width at a 
certain uj given by the slope of the band and E^', and the 
position by e{k) + T^\uj). The challenge for the analysis 
is that neither S nor e{k) is known. Nevertheless, it is 
often possible to extract these quantities by using only a 
few extra assumptions, such as particle-hole symmetry or 
the fact that T^' and T^" are related by a Kramers-Kronig 
relation. 

Here we use a particularly efficient iterative approach 
to the problem that was very recently introduced by 
Pletikosic et al. [12 . For obtaining the self-energy and 
amplitude, the data have been smoothed by downsam- 
pling the number of energy points. Figure |2] gives the 
detailed result of this method as applied to QFMLG on 
Ir(lll). The procedure determines the self-energy and 
the bare dispersion and these can be used to model the 
spectral function. The result of the model can be com- 
pared to the experimental data. Figure [2] shows this com- 
parison as well as a demonstration of the results' self- 
consistency. 

Figure [2ja) shows the MDC peak position extracted 
from the data together with the peak position obtained 
from the model and the bare dispersion. We find an excel- 
lent agreement between experimental and modelled peak 
positions. It is sufficient to model the bare dispersion by 
a second order polynomial. The obtained curvature is 
positive, as expected for the tt band of graphene below 
the Dirac point. At a binding energy of ~170 meV, the 
bare dispersion starts to deviate from the measured and 
modelled dispersion. This is the indication of a kink that 
is caused by electron-phonon coupling. The kink is also 
visible in the raw spectral function, but only barely. 

Figure |2|b) shows the energy-dependent MDC width 
obtained from the experimental data and from the model. 
The kink in the dispersion should correspond to a width 
increase at the same energy (^170 meV) but this is hard 
to discern in the data over the overall increase of the 
width in the energy range plotted here. In any event, the 
measured and modelled widths agree very well. 

According to ([T]), the self-energy also determines the 
amplitude of the spectral function and Figure [2|c) shows 
a comparison between experimental and modelled peak 
intensity. Again, the agreement is excellent. The peak 
amplitude increases notably for energies close to the 
Fermi energy and this can also be seen in the raw data. 
The peak width of the MDCs near the Fermi energy 
(~0.032 A~^) is comparable, if slightly larger, to the best 
values obtained for non-intercalated graphene on Ir(lll) 
[11120]. 

Finally, Figure |2|d) and (e) show the real and imag- 
inary part of the self-energy, respectively. Not only the 
extracted real (imaginary) part is shown, but also the 
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FIG. 2: Extraction of the self-energy S for QFMLG on 
Ir(lll).(a) Measured dispersion (MDC maximum) together 
with the model dispersion and the extracted bare dispersion. 
The inset shows a close-up of the experimental spectral func- 
tion close to the Fermi level and the spectral function calcu- 
lated from the derived bare dispersion and self-energy, (b) 
Measured MDC width together with modelled MDC width, 
(c) Measured amplitude of the spectral function together with 
modelled amplitude, (d) and (e): Demonstration of self- 
consistency, showing the modelled real part of the self-energy 
together with a Kramers-Kronig transformation (KK) of the 
imaginary part and vice versa. 



Kramers-Kronig transformed imaginary (real) part, il- 
lustrating the self-consistency of the result. The real 
part of the self-energy shows a peak at ^ 170 meV 
and so does the Kramers-Kronig transformation of the 
imaginary part. The corresponding weak step in T^" is 
only discernible in the Kramers-Kronig transformation 
of I;^ Overall, a satisfactory degree of self-consistency is 
achieved. 

The same analysis for QFMLG on SiC(OOOl) is shown 
in Fig. [3] It gives very similar results. Again, the promi- 
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FIG. 3: Extraction of the self-energy E for QFMLG on 
SiC(0001).(a) Measured dispersion (MDC maximum) to- 
gether with the model dispersion and the extracted bare 
dispersion. The inset shows a close-up of the experimental 
spectral function close to the Fermi level and the spectral 
function calculated from the derived bare dispersion and self- 
energy, (b) Measured MDC width together with modelled 
MDC width, (c) Measured amplitude of the spectral function 
together with modelled amplitude, (d) and (e): Demonstra- 
tion of self-consistency, showing the modelled real part of the 
self-energy together with a Kramers-Kronig transformation 
(KK) of the imaginary part and vice versa. 



nent structure in the data is a kink at ~ 170 meV that 
is also identified in the linewidth and in the peak inten- 
sity. The agreement between the data and the model, 
as well as the self-consistency is equally good as for the 
QFMLG on Ir(lll). The data are quite similar to the 
recent results by Forti et al [11 . 

The central result of this paper lies in the comparison 
of the results from the two QFMLG samples. Qualita- 
tively, both show the theoretically expected behaviour 
for hole-doped graphene ^ : in a region close to the 
Fermi energy, T^" is relatively constant, followed by a 



steep increase at ^ 170 meV and a linear increase at 
even higher energies. In a simple picture, this behaviour 
is explained as follows: for weakly-doped graphene, low 
energy acoustic phonons play an insignificant role for the 
electron-phonon coupling, as seen by the flat T^" near 
the Fermi level. The most significant coupling is to an 
optical phonon mode at 170 meV. The linear increase at 
even higher energies is also related to the electron-phonon 
coupling, but it stems from an increase in the electronic 
density of states away from the Dirac point, an effect that 
can have a pronounced influence on the electron-phonon 
coupling strength in semimetals [30] . 

The dominant role of the optical phonons in the scat- 
tering process is consistent with all other observations, 
such as the kink in the real part of the self-energy as well 
as the narrowing of the MDC peak and its increased in- 
tensity for binding energies smaller that 170 meV. The 
dominant coupling to optical phonons is also consistent 
with several other reported observations for electron- 
phonon coupling in graphene [5HT2]. 

Another effect that contributes to the self-energy is 
the electron-electron interaction. It can be argued, how- 
ever, that this effect is only of minor importance for 
the energy region and the kink studied here. Indeed, 
electron-electron interaction effects such as electron-hole 
and electron-plasmon processes have been shown to de- 
pend critically on the position of the Dirac point with 
respect to the Fermi level and thus on the doping |2l] . 
The observed band renormalization reported here occurs 
at a fixed binding energy of ^ 170 meV regardless of 
doping, which is seen by comparing the dispersions for 
the two substrates in Figs. [2|a) andjsj^a). Additionally, 
since the most dominant optical phonon modes are found 
at this energy in calculations 031^, we do not consider 
interactions associated with electrons to significantly in- 
fluence the observed kink. 

Note that there is an absolute offset between the dif- 
ferent MDC linewidth and T^" values for the two samples 
near the Fermi energy. The conventional interpretation 
of this would be a smaller amount of defect scattering 
for QFMLG on SiC(OOOl) than on Ir(lll). This does 
not necessarily have to be the case, however, because 
the role of not completely intercalated graphene regions 
would be quite different for the two systems. On Ir(lll), 
areas that are not completely (or not at all) intercalated 
with oxygen would still give rise to a Dirac cone, only 
at a different energy and this would broaden the spec- 
tra. On SiC, on the other hand, areas without hydrogen 
intercalation would not at all contribute to the observed 
Dirac cone. 

From the self-energy, we can now attempt a deter- 
mination of the electron-phonon mass enhancement pa- 
rameter A. As mentioned above, the key-problem with 
this type of analysis is that the determination of A be- 
comes quite uncertain in the case of weak coupling, sim- 
ply because it is difficult to quantify an effect that is 
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weak enough to be hardly observable. A far superior ap- 
proach to determine A in such a case is via a study of the 
temperature-dependent linewidth. For temperatures at 
or above the Debye temperature, this linewidth is directly 
proportional to A, permitting a straight-forward deter- 
mination of the scattering strength [151 EOl EH EH |34] . 
Unfortunately, this is not an option for QFMLG because 
the system already becomes unstable for temperatures 
far below the Debye temperature. 

We thus restrict ourselves to a simple order of magni- 
tude estimate of A, considering only the phonon mode at 
170 meV and relating the step in T," to A via [35] 
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where AH'' is the size of the jump in S'' and huph is the 
phonon energy. While it is difficult to obtain a precise 
value for AS'', an estimate consistent with both samples 
is 10 meV^ AH" < 20 meV. This leads us to estimate a 
value of A = 0.05(3). It is also possible, in principle, to 
determine A from the slope of the real part of the self- 
energy. This determination is consistent with that from 
AE" but the uncertainties are even larger. 

A value of A = 0.05(3) compares well to the few other 
studies of weakly doped graphene. Forti et al. have pub- 
lished a study for QFMLG on SiC(OOOl) and analyzed 
the MDC linewidth using a somewhat different approach 
from the one employed here. They found A ~ 0.03, 
a value consistent with ours when taking into account 
the large uncertainties of both experiments. Ulstrup 
et al have studied epitaxial graphene on Ir(lll), i.e. 
without oxygen intercalation, via temperature-dependent 
ARPES and found A ^ 9 x lO""^ [15]. 

The most important conclusion from this work is that 
quasi free-standing monolayer graphene does indeed ap- 
pear to be quasi free-standing in terms of the observed 
electron-phonon coupling: while the systems studied here 
are very different in terms of substrate materials, a wide- 
gap semiconductor vs. a metal, and intercalation ma- 
terials (H, O), the resulting electron-phonon coupling is 
similar in strength and spectral appearance. Another 
important result is that the coupling is small, consistent 
with what would be expected for pristine graphene. 
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